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Bifurcations and Strange Attractors 
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Abstract 

We reviews the theory of strange attractors and their bifurcations. All 
known strange attractors may be subdivided into the following three groups: 
hyperbolic, pseudo-hyperbolic ones and quasi-attractors. For the first ones 
the description of bifurcations that lead to the appearance of Smale- Williams 
solenoids and Anosov-type attractors is given. The definition and the descrip- 
tion of the attractors of the second class are introduced in the general case. 
It is pointed out that the main feature of the attractors of this class is that 
they contain no stable orbits. An etanol example of such pseudo-hyperbolic 
attractors is the Lorenz one. We give the conditions of their existence. In 
addition we present a new type of the spiral attractor that requires countably 
many topological invariants for the complete description of its structure. The 
common property of quasi-attractors and pseudo-hyperbolic ones is that both 
admit homoclinic tangencies of the trajectories. The difference between them 
is due to quasi-attractors may also contain a countable subset of stable peri- 
odic orbits. The quasi-attractors are the most frequently observed limit sets 
in the problems of nonlinear dynamics. However, one has to be aware that the 
complete qualitative analysis of dynamical models with homoclinic tangencies 
cannot be accomplished. 

2000 Mathematics Subject Classification: 37C29, 37C70, 37CI5, 37D45. 



1. Introduction 

One of the starling discoveries in the XX century was the discovery of dy- 
namical chaos. The finding added a fascinatingly novel type of motions — chaotic 
oscillations to the catalogue of the accustomed, in nonlinear dynamics, ones such 
as steady states, self-oscillations and modulations. Since then many problems of 
contemporary exact and engineering sciences that are modelled within the frame- 
work of differential equations have obtained the adequate mathematical description. 
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Meanwhile this also set up a question: the trajectory of what kind are to be con- 
nected to dynamical chaos in systems with 3D and higher dimensions of the phase 
space. 

The general classification of the orbits in dynamical systems is due to Poincare 
and Birkhoff. As far as chaos is concerned in particular a researcher is interested 
most of all in non-wandering trajectories which are, furthermore, to be self-liming. 
In the hierarchy of all orbits of dynamical system the high goes to the set of the 
unclosed trajectories so-called stable in Poisson's sense. The feature of them is 
that the sequence composed of the differences of the subsequent Poincare recurrence 
times of a trajectory (as it gets back into the vicinity of the initial point) may have 
no upper bound. In the case where the return times of the trajectory are bounded 
Birkhoff suggested to call such it a recurrent one. A partial subclass of recurrent 
trajectories consists of almost-periodic ones which, by definition, possess a set of 
almost-periods. Since the spectrum of the set is discrete it follows that the dynamics 
of a system with the almost-periodic trajectories must be simple. Moreover, the 
almost-periodic trajectories are also quasi-periodic, the latter ones are associated 
with the regime of modulation. The closure of a recurrent trajectory is called a 
minimal set, and that of a Poisson stable trajectory is called a quasi-minimal one. 
The last one contains also the continuum of Poisson stable trajectories which are 
dense in it. The quasi-minimal set may, in addition, contain some closed invariant 
subsets such as equilibria, periodic orbits, ergodic invariant tori with quasi-periodic 
covering and so on. All this is the reason why Poincare recurrent times of a Poisson 
stable trajectory are unbounded: it may linger in the neighborhoods of the above 
subsets for rather long times before it passes by the initial point. In virtue of such 
unpredictable behavior it seems quite reasonable the the the role of dynamical chaos 
orbits should be assigned to the Poisson stable trajectories. 

Andronov was the first who had raised the question about the correspondence 
between the classes of trajectories of dynamical systems and the types of observ- 
able motions in nonlinear dynamics in his work "Mathematical problems of auto- 
oscillations" . Because he was motivated to explain the nature of self-oscillations he 
repudiated the Poisson stable trajectory forthwith due to their irregular behavior. 
He expressed the hypothesis that a recurrent trajectory stable in Lyapunov sense 
would be almost-periodic; he also proposed to Markov to confirm it. Markov proved 
a stronger result indeed; namely that a trajectory stable both in Poisson and Lya- 
punov (uniformly) senses would be an almost-periodic one. This means that the 
Poisson stable trajectories must be unstable in Lyapunov sense to get associated to 
dynamical chaos. After that it becomes clear that all trajectories in a quasi-minimal 
set are to be of the saddle type. The importance of such quasi-minimal sets for non- 
linear dynamics as the real objects was inferred in the explicit way by Lorenz in 
1963. He presented the set of equations, known today as the Lorenz model, that 
possessed an attracting quasi-minimal set with the unstable trajectory behavior. 
Later on, such sets got named strange attractors after Rucllc and Takcns. 

Here we arrive at the following problem: how can one establish the existence 
of the Poisson stable trajectories in the phase space of a system? Furthermore, the 
applicability of the system as a nonlinear dynamics model requires that such tra- 
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jectories shall persist under small smooth perturbations. Undoubtedly the second 
problem is a way complex than that of finding periodic orbits. Below we will yield 
a list of conditions that guarantee the existence of unstable trajectories stable in 
Poisson sense. 

The most universal criterion of the existence of Poisson stable trajectories is 
the presence, in a system, of a hyperbolic saddle periodic orbit whose stable and 
unstable manifolds cross transversally along the homoclinic orbit. This structure 
implies that the set N of all of the orbits remaining in its small neighborhood 
consists of only unstable ones. Moreover, the periodic orbits are dense in N, so are 
the trajectories homoclinic to them, besides the continuum of unclosed trajectories 
stable in Poisson sense. Generally speaking wherever one can describe the behavior 
of the trajectories in terms of symbolic dynamics (using either the Bernoulli sub-shift 
or the Markov topological chains) , the proof of the existence of the Poisson stable 
trajectories becomes trivial. However, the selected hyperbolic sets by themselves 
are unstable. Nevertheless, their presence in the phase space means the complexity 
of the trajectory behavior even though they are no part of the strange attractor. 

Early sixtieth were characterized the rapid development of the theory of struc- 
tural stability initiated in the works of Anosov and Smale. Anosov was able to single 
out the class of the systems for which the hyperbolicity conditions hold in the whole 
phase space. Such flows and cascade have been named the Anosov systems. Some 
examples of the Anosov systems include geodesic flows on compact smooth man- 
ifolds of a negative curvature [3]. It is well-known that such flow is conservative 
and its set of non-wandering trajectories coincides with the phase space. Another 
example of the Anosov diffeomorphism is a mapping of an n-dimensional torus 

6 = A6 + f{6), modi, (1) 

where A is a matrix with integer entries other than 1, dct \A\ = 1, the eigenvalues 
of A do not lie on the unit circle, and f(6) is a periodic function of period 1. 

The condition of hyperbolicity of Q may be easily verified for one pure class 
of diffcomorphisms of the kind: 

6 = AO, mod 1, (2) 

which are the algebraic hyperbolic automorphisms of the torus. Automorphism J5J 
are conservative systems whose set fi of non-wandering trajectories coincides with 
the torus T" itself. 

Conditions of structural stability of high-dimensional systems was formulated 
by Smale These conditions are in the following: A system must satisfy both 
Axiom A and the strong transversality condition. 

Axiom A requires that: 

1A the non-wandering set f2 be hyperbolic; 

IB Q = Per. Here Per denotes the set of periodic points. 

Under the assumption of Axiom A the set £1 can be represented by a finite union of 
non-intersecting, closed, invariant, transitive sets fix, . . . Q p . In the case of cascades, 
any such £li can be represented by a finite number of sets having these properties 
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which are mapped to each other under the action of the diffeomorphism. The sets 
Qi, . . .Q p are called basis sets. 

The basis sets of Smale systems (satisfying the enumerated conditions) may 
be of the following three types: attractors, repellers and saddles. Repellers are the 
basis sets which becomes attractors in backward time. Saddle basis sets are such 
that may both attract and repel outside trajectories. A most studied saddle basis 
sets are one-dimensional in the case of flows and null-dimensional in the case of 
cascades. The former ones are homeomorphic to the suspension over topological 
Markov chains; the latter ones are homeomorphic to simple topological Markov 
chains [Bowen [H]]. 

Attractors of Smale systems are called hyperbolic. The trajectories passing 
sufficiently close to an attractor of a Smale system, satisfies the condition 

dist(ip(t,x),A) < fce~ At , £>0 

where k and A are some positive constants. As we have said earlier these attractors 
arc transitive. Periodic, homo- and heteroclinic trajectories as well as Poisson- 
stable ones are everywhere dense in them. In particular, we can tell one more of 
their peculiarity: the unstable manifolds of all points of such an attractor lie within 
it, i.e., W£ £ A where x G A. Hyperbolic attractors may be smooth or non-smooth 
manifolds, have a fractal structure, not locally homeomorphic to a direct product 
of a disk and a Cantor set and so on. 

Below we will discuss a few hyperbolic attractors which might be curious for 
nonlinear dynamics. The first example of such a hyperbolic attractor may be the 
Anosov torus T n with a hyperbolic structure on it. The next example of hyper- 
bolic attractor was constructed by Smale on a two-dimensional torus by means of 
a "surgery" operation over the automorphism of this torus with a hyperbolic struc- 
ture. This is the so-called DA- (derived from Anosov) diffeomorphism. Note that the 
construction of such attractors is designed as that of minimal sets known from the 
Poincarc-Donjoy theory in the case of C 1 -smooth vector fields on a two-dimensional 
torus. 

Let us consider a solid torus LI € R™, i.e., T 2 = D 2 x S 1 where D 2 is a disk 
and S 1 is a circumference. We now expand T 2 m-times (m is an integer) along 
the cyclic coordinate on S 1 and shrink it <?-times along the diameter of D 2 where 
q < 1/m. We then embed this deformated torus LIi into the original one so that its 
intersection with D 2 consists of m-smaller disks. Repeat this routine with LIi and 
so on. The set £ =G^Z 1 Hi so obtained is called a Witorius-Van Danzig solenoid. 
Its local structure may be represented as the direct product of an interval and a 
Cantor set. Smale also observed that Witorius-Van Danzig solenoids may have a 
hyperbolic structures, i.e., be hyperbolic attractors of diffcomorphisms on solid tori. 
Moreover, similar attractors can be realized as a limit of the inverse spectrum of 
the expanding cycle map |37) 

9 = m9, mod 1. 

The peculiarity of such solenoids is that they are expanding solenoids. Gener- 
ally speaking, an expanding solenoid is called a hyperbolic attractor such that its 
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dimension coincides with the dimension of the unstable manifolds of the points 
of the attractor. Expanding solenoids were studied by Williams who showed 
that they are generalized (extended) solenoids. The construction of generalized 
solenoids is similar to that of minimal sets of limit-quasi-periodic trajectories. 
Note that in the theory of sets of limit-quasi-periodic functions the Wictorius-Van 
Danzig solenoids are quasi-minimal sets. Hyperbolic solenoids arc called the Smalc- 
Williams solenoids. 

We remark also on an example of a hyperbolic attractor of a diffeomorphism 
on a two-dimensional sphere, and, consequently, on the plane, which was built by 
Plykin [23J- In fact, this is a diffeomorphism of a two-dimensional torus projected 
onto a two-dimensional sphere. Such a diffeomorphism, in the simplest case, pos- 
sesses four fixed points, moreover all of them are repelling. 

2. Birth of hyperbolic attractors 

Let us now pause to discuss the principal aspects related to the transition from 
Morse-Smale systems to systems with hyperbolic attractors. The key moment here 
is a global bifurcation of disappearance of a periodic trajectory. Let us discuss this 
bifurcation in detail following the paper of Shilnikov and Turaev |34| . 

Consider a C-smooth one-parameter family of dynamical systems in M. n . 
Suppose that the flow has a periodic orbit Lo of the saddle-node type at fi — 0. 
Choose a neighbourhood Uo of Lo which is a solid torus partitioned by the (n — 1)- 
dimensional strongly stable manifold Wff into two regions: the node region U + 
where all trajectories tend to Lq as t — > +oo, and the saddle region U~ where the 
two-dimensional unstable manifold bounded by L lies. Suppose that all of 
the trajectories of W^ return to L from the node region U + as t — * +oo and do 
not lie in W ss . Moreover, since any trajectory of W u is bi-asymptotic to Lo, W£ 
is compact. 

Observe that systems close to Xq and having a simple saddle-node periodic 
trajectories close Lq form a surface B of codimension-1 in the space of dynamical 
systems. We assume also that the family X^ is transverse to B. Thus, when fj, < 0, 
the orbit Lo is split into two periodic orbits, namely: L~ of the saddle type and 
stable L+. When \i > L disappears. 

It is clear that X^ is a Morse-Smale system in a small neighbourhood U of the 
set W u for all small /i < 0. The non-wandering set here consists of the two periodic 
orbits L+ and L~, All trajectories of U\W" tend to L+ as t — > +oo. At \i = all 
trajectories on U tend to Lo. The situation is more complex when [i > 0. 

The Poincare map to which the problem under consideration is reduced, may 
be written in the form 



where /, g and h are periodic functions of 9. Moreover, H/Hc 1 and \\h\\ci 
as fj, — > 0, m is an integer and w is a parameter defined in the set [0, 1). Diffeomor- 
phism Q is defined in a solid-torus D™ -2 x S 1 , where D"" 2 is a disk < r, r > 



x = f(x,9,fx), 

9 = m6 + g{9) + lu + h(x, 9, fi), 



mod 1, 



(3) 
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Observe that © is a strong contraction along x. Therefore, mapping |J5J is close to 
the degenerate map 

- = ' (4) 

9 = m9 + g(9)+uj, mod 1. V ' 

This implies that its dynamics is determined by the circle map 

9 = m6 + g(6) + u, modi, (5) 

where < uj < 1. Note that in the case of the flow in R 3 , the integer m may assume 
the values 0, 1. 

Theorem 2.1 If rn = and if 

max 1 1. g' (0) 1 1 < 1, 

then for sufficiently small \i > 0, the original flow has a periodic orbit both of which 
length and period tend to infinity as /j, — > 0. 

This is the "blue sky catastrophe" . In the case where m = 1 , the closure 
W£ is a two-dimensional torus. Moreover, it is smooth provided that Q is a 
diffcomorphism. In the case where m = — 1 W% is a Klein bottle, also smooth if 
is a diffcomorphism. In the case of the last theorem W£ is not a manifold. 

In the case of R" (n > 4) the constant m may be any integer. 

Theorem 2.2 Let \m\ > 2 and let \m + g'(6) > 1. Then for all fi > suffi- 
ciently small, the Poincare map 0) has a hyperbolic attractor homeomorphic to 
the Smale- Williams solenoid, while the original family has a hyperbolic attractor 
homeomorphic to a suspension over the Smale- Williams solenoid. 

The idea of the use of the saddle-node bifurcation to produce hyperbolic at- 
tractors may be extended onto that of employing the bifurcations of an invariant 
torus. We are not developing here the theory of such bifurcations but restrict ourself 
by consideration of a modelling situation. 

Consider a one-parameter family of smooth dynamical systems 

x = X(x, jLt) 

which possesses an invariant m-dimcnsional torus T m with a quasi-periodic trajec- 
tory at /i = 0. Assume that the vector field may be recast as 

V = C{n)y, 

z = fi + z 2 , (6) 
6 = 00) 

in a neighborhood of T m . Here, z G R\ y G R""™" 1 , 9 G T Tn and O(0) = 
(fix, ■ ■ ■ i fim)- The matrix C(/x) is stable, i.g., its eigenvalues lie to the left of the 
imaginary axis in the complex plane. At fi — the equation of the torus is y = 0, 
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the equation of the unstable manifold W u is y = 0, z > 0, and that of the strongly 
unstable manifold W ss partitioning the neighborhood of T™ into a node and a 
saddle region, is z — 0. We assume also that all of the trajectories of the unstable 
manifold W u of the torus come back to it as t — + +oo. Moreover they do not lie in 
W ss . On a cross-sections transverse to z = the associated Poincare map may be 
written in the form 

y = f(y,9,fj>), (7] 

6 = A6+g(6)+w + h(x,0,fj,), modi, 1 ' 

where A is an integer matrix, /, g, and h are 1-periodic functions of 9. Moreover 
ll/llc 1 ~~ *■ and \\h\\ci — *0as/x— »0, w = (u>i, ■ ■ ■ ,u) m ) where < Uk < 1. 

Observe that the restriction of the Poincare map on the invariant torus is close 
to the shortened map 

9 = A9 + g(9) + lu, modi. (8) 

This implies, in particular, that if (jHJ) is an Anosov map for all to (for example when 
the eigenvalues of the matrix A do not lie on the unit circle of the complex plane, 
and g(9) is small), then the restriction of the Poincare map is also an Anosov map 
for all /j, > 0. Hence, we arrive at the following statement 

Proposition 2.1 If the shortened map is an Anosov map for all small oj, then for 
all n > sufficiently small, the original flow possesses a hyperbolic attractor which 
is topologically conjugate to the suspension over the Anosov diffeomorphism. 

The birth of hyperbolic attractors may be proven not only in the case where the 
shortened map is a diffeomorphism. Namely, this result holds true if the shortened 
map is expanding. A map is called expanding of the length if any tangent vector 
field grows exponentially under the action of the differential of the map. An example 
is the algebraic map 

6 = A6, mod 1, 

such that the spectrum of the integer matrix A lies strictly outside the unit circle, 
and any neighboring map is also expanding. If || (G"^)) -1 1 < 1, where G = A+g(8), 
it follows then that the shortened map 

6 ==uj + A9 + g{9), modi, (9) 

is an expansion for all (i > 0. 

Shub |35| established that expanding maps are structurally stable. The study 
of expanding maps and their connection to smooth diffcomorphisms was continued 
by Williams [Williams Using the result of his work we come to the following 

result which is analogous to our theorem l2.ll namely 

Proposition 2.2 If || (G"^)) -1 1 < 1, then for all small fj, > 0, the Poincare map 
possesses a hyperbolic attractor locally homeomorphic to a direct product of K m+1 
and a Cantor set. 
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An endomorphism of a torus is called an Anosov covering if there exists a 
continuous decomposition of the tangent space into the direct sum of stable and 
unstable submanifolds just like in the case of the Anosov map (the difference is 
that the Anosov covering is not a one-to-one map, therefore, it is not a diffeomor- 
phism). The map JHJ is an Anosov covering if, we assume, | det A\ > 1 and if g(0) is 
sufficiently small. Thus, the following result is similar to the previous proposition 

Proposition 2.3 If the shortened map is an Anosov covering for all uj, then 
for all small /J, > the original Poincare map possesses a hyperbolic attractor locally 
homeomorphic to a direct product o/R m+1 and a Cantor set. 

In connection with the above discussion we can ask what other hyperbolic 
attractors may be generated from Morse-Smale systems? 

Of course there are other scenarios of the transition from a Morse-Smale sys- 
tem to a system with complex dynamics, for example, through f2-explosion, period- 
doubling cascade, etc. But these bifurcations do not lead explicitly to the appear- 
ance of hyperbolic strange attractors. 

3. Lorenz attractors 

In 1963 Lorenz ^B] suggested the model: 

x = -a{x - y), 

y = rx — y — xz, (10) 
z = —bz + xy 

in which he discovered numerically a vividly chaotic behaviour of the trajectories 
when a = 10, b = 8/3 and r = 28. the important conclusion has been derived 
from the mathematical studies of the Lorenz model: simple models of nonlinear 
dynamics may have strange attractors. 

Like hyperbolic attractors, periodic as well as homoclinic orbits are everywhere 
dense in the Lorenz attractor. Unlike hyperbolic attractors the Lorenz attractor is 
structurally unstable. This is due to the embedding of a saddle equilibrium state 
with a one-dimensional unstable manifold into the attractor. Nevertheless, under 
small smooth perturbations stable periodic orbits do not arise. Moreover, it became 
obvious that such strange attractors may be obtained through a finite number of 
bifurcations. In particular, in the Lorenz model (due to its specific feature: it has 
the symmetry group (x,y,z) *-> (— x, — y, z)) such a route consists of three steps 
only. 

Below we present a few statements concerning the description of the structure 
of the Lorenz attractor as it was done in Ej ■ The fact that we are considering 
only three-dimensional systems is not important, in principle, because the general 
case where only one characteristic value is positive for the saddle while the others 
have negative real parts, and the value least with the modulus is real, the result is 
completely similar to the three-dimensional case. Let B denote the Banach space 
of C-smooth dynamical systems (r > 1) with the C-topology, which are specified 
on a smooth three-dimensional manifold M. Suppose that in the domain U C B 
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Figure 1: Homoclinic butterfly 



each system X has an equilibrium state O of the saddle type. In this case the 
inequalities Ai < A2 < < A 3 hold for the roots A, = Xi(X), i = 1,2,3 of the 
characteristic equation at O, and the saddle value cr(X) = A2 + A3 > 0. A stable 
two-dimensional manifold of the saddle will be denoted by W s = W S (X) and the 
unstable one, consisting of O and two trajectories IY2 = Ti^(X) originating from it, 
by W u = W U (X). It is known that both W s and W u depend smoothly on X on each 
compact subset. Here it is assumed that in a certain local map V = {[x\, X2, £3)}, 
containing O, X can be written in the form 

ii = \xi + Pi(xi,X2,x 3 ), £=1,2,3. (11) 

Suppose that the following conditions are satisfied for the system Xq C U (see 
Fig.l): 

1. Ti(X ) C W S {X ), £ = 1,2 , i.e., r,(A ) is doubly asymptotic to O. 

2. Ti(Ao) and ^(Xo) approach to O tangentially to each other. 

The condition Ai < A2 implies that the non- leading manifold Wq of Wq, 
consisting of O and the two trajectories tangential to the axis x\ at the point 
O, divides Wq into two open domains: Wl and WL- Without loss of generality 
we may assume that Ti(Xo) C Wl(Xo), and hence Ti is tangent to the positive 
semiaxis xi- Let v\ and V2 be sufficiently small neighborhoods of the separatrix 
"butterfly" f = Ti U O U T2. Let A4\ stand for the connection component of the 
intersection of Wl(Xo) with Vi, which contains Ti(Xo). In the general case A4i is 
a two-dimensional C°-smooth manifold homeomorphic ether to a cylinder or to a 
Mobius band. The general condition lies in the fact that certain values Ai(Xq) and 
^2(^0)1 called the separatrix values, should not be equal to zero. 
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It follows from the above assumptions that Xq belongs to the bifurcation set 
B\ of codimension two, and B\ is the intersection of two bifurcation surfaces B\ 
and B\ each of codimension one, where B\ corresponds to the scparatrix loop 
Ti = O U r, . In such a situation it is natural to consider a two-parameter family 
of dynamical systems X(fi), fj, = (^i,yu 2 ), \lA < Mo 5 X(0) = X , such that X(fj) 
intersects with B\ only along Xq and only for /i = 0. It is also convenient to assume 
that the family X(fi) is transverse to B\. By transversality we mean that for the 
system X{p) the loop Y\(X(p)} "deviates" from W+(X(/j,)) by a value of the order 
of fii, and the loop r 2 (X(fi)) "deviates" from W+(X(fj,)) by a value of the order of 

1*2- 

It is known from |30| that the above assumptions imply that in the transition 
to a system close to Xq the separatrix loop can generate only one periodic orbit 
which is of the saddle type. Let us assume, for certainty, that the loop Ti(X ) U O 
generates a periodic orbit L\ for \i\ > and r2(Xo)UO generates the periodic orbit 
L2 for fj,2 > 0. 

The corresponding domain in U, which is the intersection of the stability 
regions for L\ and L 2 , i.e., i.e., the domain in which the periodic orbits L\ and L 2 
are structurally stable, will be denoted by Uq. A stable manifold of Li for the system 
X C Uo will be denoted by W/ and the unstable one by W". If the separatrix value 
Ai(X ) > 0, is a cylinder; if Ai(X ) < 0, W" is a Mobius band. Let us note 
that, in the case where M. is an orientable manifold, Wf will also be a cylinder 
if Ai(Xo) > 0. Otherwise it will be a Mobius band. However, in the forthcoming 
analysis the signs of the separatrix values will play an important role. Therefore, it 
is natural to distinguish the following three main cases 

Case A (orientable) Ai(X ) > 0, A 2 (X a ) > 0, 
Case B (semiorientable) Ai(X ) > 0, A 2 (X a ) < 0, 
Case C (nonorientable) Ai(X ) < 0, A 2 (X ) < 0. 

In each of the above three cases the domain Uo also contains two bifurcation surfaces 
B\ and B\: 

1. In Case A, B\ corresponds to the inclusion Y\ C W 2 and B\ corresponds to 
the inclusion T 2 C W"; 

2. In Case B, B\ corresponds to the inclusion Ti C Wf and B\ corresponds to 
the inclusion T 2 C W{; 

3. In Case C, along with the above-mentioned generated orbits L\ and L 2 , there 
also arises a saddle periodic orbit L3 which makes one revolution "along" 
Ti(Xq) and T 2 (Xq), and if both W" are Mobius bands, i = 1, 2, the unstable 
manifold W% of the periodic orbit L3 is a cylinder. In this case the inclusions 
Li C W 2 and T 2 C W§ correspond to the surfaces £?| and B\, respectively. 

Suppose that B\ and B\ intersect transversely over the bifurcational set B 2 , 
see Fig. 2. In a two-parameter family X(fi) this means that the curves B\ and B\ 
intersect at some point \i x = (fin, ^12)- Let us denote a domain lying between i3| 
and B\ by XJ\. Suppose also that for each X e U there exists a transversal D with 
the following properties: 
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Figure 2: (/ij, /i2)-bifurcation diagram 



1. The Euclidean coordinates (x, y) can be introduced on D such that 

D = {(s,y):|x|<l,|y|<2}. 

2. The equation y = describes a connection component S of the intersection 
Wq fl D such that no w-scmitrajectory that begins on S possesses any point 
of intersection with D for t > 0. 

3. The mapping T\{X) : D\ i— > Z? and Ti(X) : D 2 ^ D are defined along the 
trajectories of the system X, where 

D 1 = {(x,y):\x\<l,0<y<l}, 

D 2 = {(x iy ):\x\<l,-l>y<l}, 
and is written in the form 

x = fi(x,y), , 12 s 

where /j,^ G C r , i = 1,2. 

4. /j and <7i admit continuous extensions on S, and 

]imfi(x,y)=x**, \img l (x,y)=y**, i = 1,2. 

5. 

T X A e Pii = {(x,y) : 1/2 < x < 1, \y\ < 2}, 
T 2 D 2 e Pi 2 = : -1 < x < -1/2, |y| < 2}. 
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Let T{X) = Ti(X) | D h (/,<?) == (f l7 g t ) on D h i 
6. Let us impose the following restrictions on T(X) 

(a) ||(/*)||<1, 



1,2. 



} (13) 



(b) 1 - life)" 1 !! • H/xll > 2y/\\ (g,)- 1 1| • M • II • 

(c) IK^)- 1 !! < i, 

(d) life)- 1 - /J- H.9,11 <(l-||/ K ||)(l-||(.g,)- 1 ||). 
Hereafter, || • || = sup | • |. 

(x,y)eD\S 



It follows from the analysis of the behavior of trajectories near Wq that in a 
small neighborhood of S the following representation is valid: 



h = x\* + v>i{x,y) y a , gi 
h = x 2* + <P2{x,y)(~y) a , g 2 



yF + i/>i(x,y) y a , 

y? +Mx,y)(-y) c 



(14) 



where tpi,...,ip2 are smooth with respect to x, y for i/ / 0, and Tj(x) satisfies 
estimates l|13|l for sufficiently small y. Moreover, the limit of ipi will be denoted 
by Ai(X) and that of t/> 2 by A 2 {X). The functional Ai(X) and A 2 {X) will be 
also called the separatrix values in analogy with Ai{Xq) and ^(Xo) which were 
introduced above. Let us note that for a system lying in a small neighborhood of 
the system X all the conditions 1-6 are satisfied near S. Moreover, the concept of 
orientable, semiorientable and nonorientable cases can be extended to any system 
X G U. It is convenient to assume, for simplicity, that A 12 {X) do not vanish. It 
should be also noted that the point Pj with the coordinates (x**,y**) is the first 
point of intersection of Ti(X) with D. 
Let us consider the constant 



q = 



i h- ll^HIICgy)- 1 !! -t- a/i — llCg^)- 1 !! 2 !!^)!! — ^ik.^)- 1 !! llg^ll II (g^)- 1 ^ 

silCg,)- 1 !! 



(15) 



Conditions (|15|l implies that q > 1 and hence all the periodic points will be of the 
saddle type. 

Let E denote the closure of the set of points of all the trajectories of the 
mapping T(X), which are contained entirely in D. S is described most simply in 
the domain U\. Here the following theorems hold. 

Theorem 3.1 If X G U\, T(X) \ £ is topologically conjugated with the Bernoulli 
scheme (a, Q 2 ) with two symbols. 



Theorem 3.2 The system X G U 2 , has a two-dimensional limiting set Q, which 
satisfies the following conditions: 
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1. Q is structurally unstable. 

2. [r 1 ur 2 no]ca 

3. Structurally stable periodic orbits are everywhere dense in fi. 

4- Under perturbations of X periodic orbits in Q disappear as a result of matching 
to the saddle separatrix loops T\ and 

Note that in this case the basic periodic orbits will not belong to il. In terms 
of mappings, the properties of f2 can be formulated in more detail. Let us first 
single out a domain D on D as follows: we assume that in Case A 

D ={(x,y) G L>i U L> 2 | y 2 (x) <y< 

and in Case B 

D ={(x,y) G D x UD 2 | y 12 {x) <y< 

where y = 2/12(2;), \x\ < 1, denotes a curve in D whose image lies on the curve 
y = yi(x); and finally in Case C 

The closure of points of all the trajectories of the mapping T(X), which are 
entirely contained in D, we will denote by E. 

Theorem 3.3 Let X G U 2 - Then: 

I. £ is compact, one- dimensional and consists of two connection components in 
Cases A and C, and of a finite number of connection components in Case B. 
II. D is foliated by a continuous stable foliation H + into leaves, satisfying the 
Lipschitz conditions, along which a point is attracted to E; inverse images of 
the discontinuity line S : y = (with respect to the mapping T k , k = 1, 2, . . .) 
are everywhere dense in D. 

III. There exits a sequence of T(X) -invariant null- dimensional sets Ak, k G Z + , 
such thatT(X) \ is topologically conjugated with a finite topological Markov 
chain with a nonzero entropy, the condition G A^ +1 being satisfied, and 
Afe — > E as k — > oo. 

IV. The non-wandering set E x e E is a closure of saddle periodic points of T(X) 
and either Ei = E or Si = E + U E~, where: 

1. E~ is null- dimensional and is an image of the space Q~ of a certain TMC 
(G~, a) under the homeomorphism [3 : £~ i— ► E~ which conjugates a \ 
and T(X) | 

l(X) 

£"= U £-,Z(X)<oo, 

where 

T(X)£- = E" , E" 1 n E" 2 = 

for m\ 7^ ni2 and T(X) \ E~ is transitive; 

2. E + is compact, one- dimensional and 

3. if E + n E~ = 0, E + is an attracting set in a certain neighbourhood; 
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4- j/S + HS ^ 0, thenYj + C\Tj = S+nE m for a certain m, and this intersection 
consists of periodic points of no more than two periodic orbits, and 

(a) if^m is finite, E + is lu -limiting for all the trajectories in a certain neighbour- 
hood; 

(b) if E~ is infinite, E + is not locally maximal, but is uj-limiting for all the 
trajectories in D, excluding those asymptotic to E _ \S + . 

Below we will give the conditions under which the existence of the Lorenz 
attractor is guaranteed. 

Consider a finite- number parameter family of vector field defined by the system 
of differential equations 

x = X[x, fx), (16) 

where x £ R n+1 , /i £ R m , and X{x, fx) is a C r -smooth functions of x and /Lt. Assume 
that following two conditions hold 

A. System <|16|) has a equilibrium state O(0, 0) of the saddle type. The eigenvalues 
of the Jacobian at 0(0,0) satisfy 

ReA n < ■ • • RcA 2 < Ai < < A . 

B. The separatrices IT and I2 of the saddle O(0,0) returns to the origin as 
t — > +00. 

Then, for /i > in the parameter space there exists an open set V, whose 
boundary contains the origin, such that in V system (I16H possesses the Lorenz 
attractor in the following three cases |31| : 
Case 1. 

A Ti and T return to the origin tangcntially to each other along the dominant 
direction corresponding to the eigenvalue Ai; 

B 

1 Ai ReXi 

- < 7 < 1, Vi > 1, 7 = - — , Vi = — ; 

A The separatrix values A\ and A2 (see above) are equal to zero. 

In the general case, the dimension of the parameter space is four since we may 
choose /ii^j to control the behaviour of the separatrices Ti^ and /X3.4 = A3. 4. In 
the case of the Lorenz symmetry, we need two parameters only. 
Case 2. 

A Ti and T belong to the non-leading manifold W ss £ W s and enter the saddle 
along the eigen-direction corresponding to the real eigenvector A2 

B 

1 Ai RcAj 

- < 7 < 1, vi > 1, 7 = - — , Vi = — ; 

l Ao Ao 

In the general case, the dimension of the phase space is equal to four. Here, ^3,4 
control the distance between the separatrices. 
Case 3. 
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a r 1)2 i w ss - 

B 7 = 1; 

C A ia ^ 0, and \A ia \ < 2. 

In this case m = 3, ^3 = 7 — 1. 

In the case where the system is symmetric, all of these bifurcations are of 
codimension 2. 

In A. Shilnikov |57||2S] it was shown that both subclasses (A) and (C) are real- 
ized in the Shimizu-Morioka model in which the appearance of the Lorenz attractor 
and its disappearance through bifurcations of lacunae are explained. Some systems 
of type (A) were studied by Rychlik [23] and those of type (C) by Robinson |24| . 

The distinguishing features of strange attractors of the Lorenz type is that 
they have a complete topological invariant. Geometrically, we can state that two 
Lorenz-like attractors are topologically equivalent if the unstable manifolds of both 
saddles behave similarly. The formalization of "similarity" may be given in terms 
of kneading invariants which were introduced by Milnor and Thurston while 
studying continuous, monotonic mappings on an interval. This approach may be 
applied to certain discontinuous mappings as well. Since there is a foliation (see 
above) we may reduce the Poincare map to the form 

x = F(x,y), , , 

V = G(y), [U) 

where the right-hand side is, in general, continuous, apart from the discontinuity 
line y = 0, and G is piece- wise monotonic. Therefore, it is natural to reduce l|17ll 
to a one-dimensional map 

V = G(y), 

by using the technique of taking the inverse spectrum, Guckcnhcimcr and Williams 
|16j showed that a pair of the kneading invariants is a complete topological invariant 
for the associated two-dimensional maps provided inf |G"| > 1. 



4. Wild strange attractor 

In this section, following the paper by Shilnikov and Turaev we will 

distinguish a class of dynamical systems with strange attractors of a new type. 
The peculiarity of such an attractor is that it may contain a wild hyperbolic set. 
We remark that such an attractor is to be understood as an almost stable, chain- 
transitive closed set. 

Let X be a smooth (C r , r > 4) flow in R n (n > 4) having an equilibrium state 
O of a saddle-focus type with characteristic exponents 7, — A± iuj, —a\, • • • , — a„_3 
where 7 > 0, < A < Rca.,, w^fl. Suppose 

7 > 2A. (18) 

This condition was introduced in |21| where it was shown, in particular, that it 
is necessary in order that no stable periodic orbit could appear when one of the 
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separatrices of O returns to O as t — > +00 (i.e., when there is a homoclinic loop; 
see also [221 ) - 

Let us introduce coordinates (x, y, z) (x € R 1 , y £ R 2 , z £ R n ~ 3 ) such that 
the equilibrium state is in the origin, the one-dimensional unstable manifold of O 
is tangent to the x-axis and the (n — l)-dimensional stable manifold is tangent 
to {x = 0}. We also suppose that the coordinates correspond to the leading 
exponents A± iu> and the coordinates z correspond to the non-leading exponents a. 

Suppose that the flow possesses a cross-section, say, the surface II : {|| y || = 
1,|| z ||< 1,| x |< 1}. The stable manifold W s is tangent to {x = 0} at O, 
therefore it is locally given by an equation of the form x = h s (y,z) where h s is a 
smooth function h s (0, 0) = 0, (h s )'(0,0) = 0. We assume that it can be written 
in such form at least when (|| y ||< 1, j z j|< 1) and that | h s |< 1 here. Thus, 
the surface II is a cross-section for W* oc and the intersection of Wf oc with II has 
the form IIq : x = ho((p,z) where tp is the angular coordinate: y\ =\\ y \\ cos</j, 
2/2 =|| y || simp, and fig is a smooth function — 1 < ho < 1. One can make ho = 
by a coordinate transformation and we assume that it is done. 

We suppose that all the orbits starting on n\IIo return to II, thereby defining 
the Poincarc map: T+ : 11+ -> II T_ : n_ -»■ II, where 11+ = II n {x > 0} 
II_ =nn{i<0}. It is evident that if P is a point on II with coordinates [x, (p, z), 
then 

lim T.{P)=P 1 _, lim T+ (P) = Pj , 

where P 1 — and P+ are the first intersection points of the one-dimensional separa- 
trices of O with II. We may therefore define the maps T+ and T_ so that 

r_(n ) = Pl, r+(n ) = p\. (19) 

Evidently, the region V filled by the orbits starting on II (plus the point O 
and its separatrices) is an absorbing domain for the system X in the sense that the 
orbits starting in dT> enter V and stay there for all positive values of time t. By 
construction, the region V is the cylinder {|| y \\< 1, || z \\< 1, | x |< 1} with two 
glued handles surrounding the separatrices, see Fig. 3. 

We suppose that the (semi)flow is pseudohyperbolic in T>. It is convenient for 
us to give this notion a sense more strong than it is usually done Namely, we 
propose the following 

Definition. A semi-flow is called pseudohyperbolic if the following two conditions 
hold: 



A At each point of the phase space, the tangent space is uniquely decomposed 
(and this decomposition is invariant with respect to the linearized semi-flow) 
into a direct sum of two subspaces N± and N2 (continuously depending on 
the point) such that the maximal Lyapunov exponent in Ni is strictly less 
than the minimal Lyapunov exponent in N2: at each point M , for any vectors 
u G Ni(M) and v G N 2 (M) 

,. i, Kll . f 1, IMI 

lim sup — In — — — < lim mi — in ■ 



t^+oc t \\u\\ i^ + oc: t \\v\ 
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Figure 3: Construction of the wild attractor 

where u t and v t denote the shift of the vectors u and v by the semi-flow lin- 
earized along the orbit of the point M ; 

B The linearized flow restricted on N2 is volume expanding: 

V t > const ■ e at V 

with some a > 0; here, Vq is the volume of any region in N2 and Vt is the 
volume of the shift of this region by the linearized semi-flow. 

The additional condition B is new here and it prevents of appearance of sta- 
ble periodic orbits. Generally, our definition includes the case where the maximal 
Lyapunov exponent in Ni is non-negative everywhere. In that case, according to 
condition A, the linearized semi- flow is expanding in 7V2 and condition B is satisfied 
trivially. In the present paper we consider the opposite case where the linearized 
semi-flow is exponentially contacting in Ni , so condition B is essential here. 

Note that the property of pseudo-hypcrbolicity is stable with respect to small 
smooth perturbation of the system: according to |17j the invariant decomposition 
of the tangent space is not destroyed by small perturbations and the spaces Ni and 
N2 depend continuously on the system. Hereat, the property of volume expansion 
in N2 is also stable with respect to small perturbations. 

Our definition is quite broad; it embraces, in particular, hyperbolic flows for 
which one may assume (N!,N 2 ) = {N s , N u © N ) or (N U N 2 ) = (N s © N ,N U ) 
where N s and N u are, respectively, the stable and unstable invariant subspaces and 
No is a onc-dimcnsional invariant subspacc spanned by the phase velocity vector. 
The geometrical Lorenz model from |21 EJ or ^Sj belongs also to this class: here 
Ni is tangent to the contracting invariant foliation of codimension two and the 
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expansion of areas in a two-dimensional subspace N 2 is provided by the property 
that the Poincare map is expanding in a direction transverse to the contracting 
foliation. 

In the present paper we assume that N\ has codimension three: dimN\ = n — 3 
and dimN2 = 3 and that the linearized flow (at t > 0) is exponentially contracting 
on Nx- Condition A means here that if for vectors of N 2 there is a contraction, it 
has to be weaker than those on Ni. To stress the last statement, we will call Nx 
the strong stable subspace and N 2 the center subspace and will denote them as N ss 
and N c respectively. 

We also assume that the coordinates (x, y, z) in R n are such that at each point 
of T> the space N ss has a non-zero projection onto the coordinate space z, and N c 
has a non-zero projection onto the coordinate space (x,y). 

Note that our pscudohypcrbolicity conditions are satisfied at the point O from 
the very beginning: the space N ss coincides here with the coordinate space z, 
and N c coincides with the space (x,y); it is condition l|18fl which guarantees the 
expansion of volumes in the invariant subspace (x,y). The pseudohyperbolicity of 
the linearized flow is automatically inherited by the orbits in a small neighborhood 
of O. Actually, we require that this property would extend into the non-small 
neighborhood T> of O. 

According to the exponential contraction in N ss implies the existence of 
an invariant contracting foliation J\f ss with C r -smooth leaves which are tangent to 
N ss . As in 0, one can show that the foliation is absolutely continuous. After a 
factorization along the leaves, the region T> becomes a branched manifold (since 
T> is bounded and the quotient-semiflow expands volumes it follows evidently that 
the orbits of the quotient-semiflow must be glued on some surfaces in order to be 
bounded; cf . |39j ) . 

The property of pscudohypcrbolicity is naturally inherited by the Poincare 
map T = (T+, TJ) on the cross-section II: here, we have: 

A* There exists a foliation with smooth leaves of the form (x, tp) = h(z) |-i< z <i, 
where the derivative h'(z) is uniformly bounded, which possesses the following 
properties: the foliation is invariant in the sense that if I is a leaf, then 
TI 1 (I n T+ (11+ U n ) ) and TZ 1 (I n T_ (n_ U n ) ) are also leaves of the foliation 
(if they are not empty sets); the foliation is absolutely continuous in the sense 
that the projection along the leaves from one two-dimensional transversal to 
another increases or decreases the area in a finite number of times and the 
coefficients of expansion or contraction of areas are bounded away from zero 
and infinity; the foliation is contracting in the sense that if two points belong 
to one leaf, then the distance between the iterations of the points with the 
map T tends to zero exponentially; 

B* The quotient maps T + and T_ arc area-expanding. 

Statement 3.1 Let us write the map T as 



(x,<p) = g(x,(p,z), z = f(x,<p,z), 
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where f and g are functions smooth at x ^ and discontinuous at x = 0: 



Let 



lim n (ff,/) = {x-,tp-,z-) = P±, Urn (g, f) = (x+,tp + ,z+) = P+ . 

x—*— x— >+0 



dct ± . (20) 

9(x,<p) 



Denote 



If 



\ = df df ( dg \ 1 dg g = gj f Og 

C — ( dg \ dg D — ( d 9 \ 



limC= 0, lim || A II II D 11= 0, (21) 

X— »0 X— >Q 



sup VII ^11 11^11 + J sup || P || sup || C || < 1, (22) 
Pen\n y Pen\n Pen\n 

then the map has a continuous invariant foliation with smooth leaves of the form 
(x,tp) — h(z) |-i< 2 <i where the derivative h'(z) is uniformly bounded. If, addition- 
ally, 

sup || A || + / sup || P |j sup || C || < 1, (23) 
pen\n y Pen\n Pen\n 

then the foliation is contracting and if, moreover, for some f3 > 



the functions A | x | D \ x \^ , P, C are uniformly bounded and 
Holder continuous, 

and and — — — D \ x \@ are uniformly bounded, 

oz o(x, ip) 

then the foliation is absolutely continuous. The additional condition 



(24) 



sup Vdct D + / sup || P || sup || C || < 1 (25) 
Pen\n y Pen\n Pen\n 

guarantees that the quotient map T expands areas. 

It follows from |^|22j that in the case where the equilibrium state is a saddle- 
focus, the Poincare map near IIo = II n W s is written in the following form under 
some appropriate choice of the coordinates. 

(x,0) = Q±(Y,Z), z = R±{Y,Z). (26) 

Here 



Y-\t\p ( co < nin \ x \ +l P) sin(Qln \x\+ip) 

11 ^ -8in(nin|x|+v>) cos(filn | x \ +tp) ' + U j ' 

Z = $ 2 (x,(p,z), 



(27) 
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where p = A/7 < 1/2 (see JTHJ0) ^ = w /7 an d) for some n > p, 

11^^11=0(1,1-), 0<p+U|<r-2; (28) 

the functions Q±, i?± in 126|) ("+" corresponds to x > - the map T+, "-" corre- 
sponds to x < - the map T_) are smooth functions in a neighborhood of (Y, Z) = 
for which the Taylor expansion can be written down 

Q± = 0±, <p±) + a±Y + b±Z + ■ ■ ■ , R± = z± + c ± Y + d±Z + ■■■ . (29) 

It is seen from (|26[) - l(29[) that if O is a saddle-focus satisfying (|18J) . then if 
a + 7^ and a_ 7^ 0, the map T satisfies conditions l|21|l and (|24|> with j3 £ [p,r])- 
Furthermore, analogues of conditions l|20|) . i|22[l . (|23|) . (|25|l arc fulfilled where the 
suprcmum should be taken not over | x |< 1 but it is taken over small x. The 
map l(26|) . (|27(l . 1)29(1 is easily continued onto the whole cross-section II so that the 
conditions of the lemma were fulfilled completely An example is given by the map 

x = 0.9 I x \ p cos(ln I x | +tp), 

<p = 3\x\P sin(ln | x \ +<p), (30) 
z = (0.5 + 0.12 I x Hsignx 

where 0.4 = p < rj 

As stated above, the expansion of volumes by the quotient-semifiow restricts 
the possible types of limit behavior of orbits. Thus, for instance, in T> there may be 
no stable periodic orbits. Moreover, any orbit in T> has a positive maximal Lyapunov 
exponent. Therefore, one must speak about a strange attractor in this case. 

Beforehand, we recall some definitions and simple facts from topological dy- 
namics. Let X t P be the time-i shift of a point P by the flow X. For given e > 
and r > let us define as an (e, r)-orbit as a sequence of points Pi, P2, ■ ■ ■ ,Pk such 
that Pi+i is at a distance less than e from X t Pi for some t > t. A point Q will be 
called (e,t)- attainable from P if there exists an (e,r)-orbit connecting P and Q; 
and it will be called attainable from P if, for some r > 0, it is (e, r)-attainable from 
P for any e (this definition, obviously, does not depend on the choice of r > 0). A 
set C is attainable from P if it contains a point attainable from P. A point P is 
called chain- recurrent if it is attainable from X t P for any t. A compact invariant 
set C is called chain-transitive if for any points P € C and Q £ CC and for any 
e > and r > the set C contains an (e, r)-orbit connecting P and Q. Clearly, all 
points of a chain-transitive set are chain-recurrent. 

A compact invariant set C is called orbitally stable, if for any its neighborhood 
U there is a neighborhood V(C) C [7 such that the orbits starting in V stay in 
[7 for all t > 0. An orbitally stable set will be called completely stable if for any 
its neighborhood J7(C) there exist £0 > 0, r > and a neighborhood V(C) C [7 

such that the (eoj r)-orbits starting in V never leave [/. It is known, that a set C 

00 

is orbitally stable if and only if C = O Uj where {Uj}°^L 1 is a system of embedded 
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open invariant (with respect to the forward flow) sets, and C is completely stable 
if the sets Uj are not just invariant but they are absorbing domains (i.e.; the orbits 
starting on dUj enter inside Uj for a time interval not greater than some Tf, it 
is clear in this situation that (e, T)-orbits starting on dUj lie always inside Uj if 
e is sufficiently small and r > Tj). Since the maximal invariant set (the maximal 
attractor) which lies in any absorbing domain is. evidently asymptotically stable, 
it follows that any completely stable set is cither asymptotically stable or is an 
intersection of a countable number of embedded closed invariant asymptotically 
stable sets. 

Definition. We call the set A of the points attainable from the equilibrium state O 
the attractor of the system X . 

This definition is justified by the following theorem. 

Theorem 3.4 The set A is chain-transitive, completely stable and attainable from 
any point of the absorbing domain T>. 

Let us consider a one-parameter family of such systems assuming that: 
a homoclinic loop of the saddle-focus O exists at /i = 0, 

i.e., one of the separatrices (say, r + ) returns to O as t — > +oo. In other words, 
we assume that the family X^ intersects, at ju = 0, a bifurcational surface filled 
by systems with a homoclinic loop of the saddle-focus and we suppose that this 
intersection is transverse. The transversality means that when fj, varies, the loop 
splits and if M is the number of the last point of intersection of the separatrix T + 
with the cross-section II at \i = (P^ G n at (J, = 0), then the distance between 
the point Pjjjj- and IIo changes with a "non-zero velocity" when [i varies. We choose 
the sign of \i so that P^ G 11+ when fx > (respectively, P^ { G II _ when fi < 0). 

Theorem 3.5 There exists a sequence of intervals Aj (accumulated at \i = 0) such 
that when /i G Aj, the attractor A^ contains a wild set (non-trivial transitive closed 
hyperbolic invariant set whose unstable manifold has points of tangency with its 
stable manifold). Furthermore, for any fi* G Aj, for any system C '-close to a 
system X* , its attractor A also contains the wild set. 

We have mentioned earlier that the presence of structurally unstable (non- 
transverse) homoclinic trajectories leads to non-trivial dynamics. Using results 
[111 IT3"] we can conclude that the systems whose attractors contain structurally 
non-transverse homoclinic trajectories as well as structurally unstable periodic or- 
bits of higher orders of degeneracies are dense in the given regions in the space of 
dynamical systems. In particular, the values of n are dense in the intervals Aj for 
which an attractor of the system contain a periodic orbit of the saddle-saddle type 
along with its three-dimensional unstable manifold. For these parameter values, the 
topological dimension of such an attractor is already not less than three. The latter 
implies that the given class of systems is an example of hyperchaos. 

5. Summary 

The above listed attractors arc, in the ideal, the most suitable images of dy- 
namical chaos. Even though some of them arc structurally unstable, nevertheless it 
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is important that no stable periodic orbits appear in the system under small smooth 
perturbations. Nonetheless, excluding the attractors of the Lorenz type, no others 
have been ever observed in nonlinear dynamics so far. A research has frequently 
to deal with the models in which despite the complex behavior of the trajectories 
appears to be so visually convincing, nevertheless explicit statements regarding the 
exponential instability of the trajectories in the limit set can be debatable, and 
therefore should be made with caution. In numeric experiments with such model 
one finds a positive Lyapunov exponent, a continuous frequency spectrum, fast de- 
caying correlation functions etc., i.e. all the attributes of dynamical chaos seem to 
get fulfilled so that the presence of the dynamical chaos causes no doubts. However, 
this "chaotic attractor" may and often do contain countably many stable periodic 
orbits which have long periods and week and narrow attraction basins. Besides 
the corresponding stability regions are relatively miniature in the parameter space 
under consideration and whence those orbits do not reveal themselves ordinarily 
in numeric simulations except some quite large stability windows where they are 
clearly visible. If it is the case, the quasi-attractor 0] is a more appropriate term for 
such chaotic set. The natural cause for this rather complex dynamics is homoclinic 
tangencics. Today the the systems with homoclinic tangencies are the target of 
many studies. We briefly outline some valuable facts proven for 3D systems and 
2D diffcomorphisms. It will be clear that these results will also hold for the general 
case where there may be some other peculiarities, as for instance the co-existence 
of countable sets of saddle periodic orbits of distinct topological types (see 

We suppose that the system possesses an absorbing area embracing the hyper- 
bolic basis set in which the stable and unstable subsets may touch each other. If 
it is so, such a hyperbolic set is called wild. It follows then that either the system 
itself or a close one will have a saddle periodic orbit with non-transverse homoclinic 
trajectory along which the stable and unstable manifolds of the cycle have the tan- 
gency. In general, the tangency is quadratic. Let the saddle value | A 'y | be less then 
1, where A and 7 are the multipliers of the saddle periodic orbit. This condition 
is always true when the divergence of the vector field is negative in the absorbing 
area. Therefore, near the given system there will exist the so-called Newhouse re- 
gions (20] in the space of the dynamical system, i.e. the regions of dense structural 
instability. Moreover, a system in the Newhouse region has countably many stable 
periodic orbits which cannot principally be separated from the hyperbolic subset. 
If additionally this hyperbolic set contains a saddle periodic orbit with the saddle 
value exceeding one, then there will be a countable set of repelling periodic orbits 
next to it, and whose closure is not separable from the hyperbolic set either. The 
pictures becomes ever more complex if the divergence of the vector field is sign- 
alternating in the absorbing area. Such exotic dynamics requires infinitely many 
continuous topological invariant — moduli, needed for the proper description of the 
system in the Newhouse regions. This result comes from the fact that the systems 
with the countable set of periodic orbits of arbitrary high degrees of degeneracies 
are dense in the Newhouse regions ^2 EH] • That is why we ought to conclude in a 
bitter way: the complete theoretical analysis of the models, which admit homoclinic 
tangencies, including complete bifurcation diagrams and so forth is non realistic. 
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